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A relativistic quark potential model is used to do a systematic search for quasi-stable dibaryon 
states in the u, d, and s three flavor world. Flavor symmetry breaking and channel coupling effects 
are included and an adiabatic method and fractional parentage expansion technique are used in the 
calculations. The relativistic model predicts dibaryon candidates completely consistent with the 
nonrelativistic model. 
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I. INTRODUCTION 



Since Jaffe predicted the first dibaryon, the H particle m, there have been many experimental and theoretical 
efforts to search for dibaryons. Up to now, there are no experimentally well established dibaryon states except the 
deuteron. Theoretically almost all QCD models, including lattice QCD, predict that there should be dibaryon states 
M . However no model has yet achieved an acceptable level of quantitative reliability. Either the model Hamiltonian 
is oversimplified or the model Hilbcrt space is rather restricted, or both. We developed a nonrelativistic model, which 
we termed the quark derealization color screening model(QDCSM) J3|-||]. The model Hamiltonian and Hilbert space 
are reasonable enough to yield a qualitatively correct N-N interaction including both repulsive core and intermediate 
range attraction and to fit all of the four (JJ)=(0f) (fO) (00) and (11) channel phase shifts qualitatively. This model 
reproduces an almost correct deuteron state and verifies dynamically that there are two kinds of dibaryons One 
kind is a loosely bound, low spin (J), two baryon state with only slight quark derealization, of which the deuteron is 
prototypical. Another kind is a tightly bound, high spin, six quark state with significant quark derealization; here, 
the d*(IJ=03) [|| is a prime example. (We do not discuss, however, potential dibaryons with more involved internal 
structure, such as those with non- valence internal structure, such as would correspond to NN-7T bound states [Q, or 
with orbital excitations of the valence quarks §.) 

It is believed that systems of light (u,d and s) quarks require relativistic dynamics. Therefore it is better to have 
a relativistic model calculation to check if those dibaryon states predicted by the nonrelativistic model(QDCSM) are 
robust against relativistic effects. A relativistic quark potential model(LAMP) || has been developed for this purpose 
and the fractional parentage technique developed for the nonrelativistic dibaryon calculation has been extended to 
the relativistic case as well jlO]]. Together, these make a systematic relativistic quark potential model dibaryon search 
feasible and the results are reported here. 



We neglect the small current quark mass for u and d and treat them as massless, m u — = m = 0. For a single 
baryon, we assume a relativistic Hamiltonian with a scalar linear confinement ||, 



II. THE RELATIVISTIC QUARK POTENTIAL MODEL 
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(1) 



Hi = di-pi+ f3i(mi + V(ri)) 



(2) 



V( n ) = k 2 



(n - r ) 



(3) 



Hij = g s {mi)g s {mj) 




(4) 
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where a,,/?, are Dirac matrices, Xi is the color SU(3) generator, rrij is the quark mass. V(r) is a phenomenological 
confinement potential, which we assume is a Lorentz scalar and r,; is the modulus of the three-coordinate of the z-th 
quark relative to the origin. The color Coulomb interaction due to gluon exchange is assumed to be absorbed in the 
form of Eq.(3) even though it is a Lorentz vector. This is done to simplify the numerical calculation. The parameters 
fc 2 =0.9 GeV/fm and tq = 0.57 fm are fixed by the condition that the eigenenergy E of the following Dirac equation 
& 

Hi^in) = Eip^in) (5) 

is one third of the average of the nucleon (N) and Delta(A) masses, E = ^(N + A), for massless quarks. The lowest 
energy wavefunction solution is: 

Mr) = P u {kr)e~ k2r2 l 2 

Mr)=Pi(kr)e- k2r2 / 2 (6) 
P u (kr) = 1 + 0.04287fcr + 0.00457(fcr) 2 + • • • 

Pi(kr) = -fcr(0.46330 - 0.08767fcr + • • •) 

Hij is the color magnetic interaction due to gluon exchange, assumed to keep the form of single gluon exchange. A 
Gaussian gluon propagator is adopted for the color magnetic interaction to simulate the confinement property of a 
color gluon on the one hand and to simplify the numerical calculation on the other hand, v — ( °' 2 ^. eV ) 2 = 1.0 fm~ 2 
is chosen to be about the size of a hadron. The combination g 2 (m — Q)A — —1.350 GeV is fixed by N — A mass 
difference, where A is the matrix element of the spatial part of the (color) current-current interaction matrix element. 
Note that the quark-gluon effective coupling is assumed to vary with a scale related to the quark mass. 

Note that, having so-modified the form of single gluon exchange, this two-body-interaction Hamiltonian could 
equally well be taken to represent the effects of instantons pl| . The color structure, limited range, and quark-mass- 
dependent effective strength are all also features found in the effective quark-quark interactions produced by the 
propagation of quarks in the presence of instantons. The diquark correlations expected in instanton models are thus 
implicitly implemented by this ify. 

The single baryon ground state wave function (WF) is assumed to be a product of single quark WFs, 
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^ B (123) = X c(123) 



(7) 



where x c (123) is a three quark color singlet state. r\f i is the single quark flavor WF. [ \su means the individual quark 
spin-flavors are coupled to total baryon strangness S, isospin / and spin J. 

The strange quark mass m s — 307 MeV is determined by the A — TV mass difference. (Although a factor of about 
two larger than conventional in quark models, this difference can be traced to a virial theorem factor of one-half which 
arises from the use of a Dirac scalar potential.) The ratio, g s (m s )/g s (Q), is determined by the overall fit of the octet 
and decuplet baryon masses. The fitted baryon masses and the value of g s (m s )/g s (0) are shown in Table 1. 

Table 1. Baryon masses. g s (m s )/g s (0)=l/l.26 



N A A £ 3 E* E* fl 
exp. 939 1232 1116 1193 1318 1385 1533 1675 
theor. 939 1232 1116 1182 1333 1376 1527 1683 



For a two baryon system, we change the Hamiltonian slightly from that given in Eq.(l): i now runs from 1 to 6 
and the confinement potential is replaced by 



(8) 



where R\ and R 2 are two baryon 'centers' || and V(i~i) is given in Eq.(3). This confinement potential is devised to 
simulate the effect that, as the two baryons are brought close together, the quark matter density in between them will 
increase gradually causing the QCD vacuum to change (also gradually) from nonperturbative to perturbative. As 
one quark is removed to a distance large compared to the separation of the two centers, the confining potential should 
be the same as for an isolated baryon. These effects are achieved by truncating the value of the confining potential 
on the midplane between the two baryon centers, R\ and R 2 . 

The two baryon WF is assumed to be a linear combination of the following Dirac cluster WFs, 



(9) 



^(Bi) = Xci (123) 
V(B 2 ) = Xc 2 (456) 



-I S1/1J1 



n<^)%(j) 

J'=4 



(10) 



S2I2J2 



^{n) = (^M - £1) + e{x)ip ai (fi - R 2 ))/N(x) 
< (fj) = (VV, ft - R2) + e(x)ip aj (fj - Ri))/N(x) 



(11) 



x = \R 1 -R 2 \ 
N 2 (x) = 1 + € 2 {x) + 2e(x) (i/) at (fi - RiMaAn - R 2 ) 



(12) 

A is the normalized antisymmetrization operator, a = (SIJ) and K represent the other quantum numbers related 
to B\ and B 2 , Xcn and % C2 are both the three quark color singlet states. 



A physical two baryon state with quantum number a = (SIJ) is 

f a (9 6 ) = ^C aJf * aJ f(l 6 ) (13) 

K 

The summation K runs over all the possible two baryon states consisting of the octet and decuplet baryons listed in 
Table I. This is sufficient as it spans the space of relevant asymptotic states. 

The six quark Hamiltonian is first diagonalized in the Dirac cluster basis space given by Eq.(f3). To simplify the six 
quark Hamiltonian matrix element calculations, a relativistic extension of the fractional parentage expansion method 
is used [ [l0| . The one-body and two-body matrix elements are calculated by a program developed in the preliminary 
relativistic dibaryon calculations 

The lowest eigen-energy obtained in the diagonalization is x and e(x) dependent, and e(x) is varied at each x 
to obtain the minimum. The difference between the minimum at x and infinite separation is assumed to be the 
effective potential of the two baryons. (The minimum at x = 3 fm is already negligibly different from values at 
larger separation.) We expect that spurious effects such as center-of-mass (cm.) motion largely cancel out in the 
effective potential due to the process described above. Therefore, we have neither subtracted the cm. motion energy 
separately for a single baryon nor for the two baryons together. Following Ref. ||, we use the following expression to 
calculate the dibaryon masses, 

o f- 2 

M a {q 6 ) = (Mi + M 2 ) a + mm[V a {x)\ + j (14) 

where min[V^.(x)] is the minimum of the effective interaction V a (x), and (Mi + M2) a is the channel weighted experi- 
mental two baryon mass, 

(Mi + M 2 ) a = \C aK \ 2 < y aK (q 6 )\y a K(q 6 ) > (Mi + M 2 ) aK , (15) 

K 

at the point xq at which the minimum occurs. 

A zero point oscillation energy of ^J 1 x -i has been taken into account as explained in Ref. ||. Even though the 
internal motion of the quarks is relativistic, we assume the relative motion of the two baryon centers still can be 
approximated as a nonrelativistic oscillation around the equilibrium separation xq , where /i a is the weighted reduced 
channel mass, 

Ma - £ \C aK \ 2 < * aK (q 6 )\* a K(q b ) > ( ) • (16) 

K \M 1 +M 2 J aK 

This model has been used to calculate the binding energies of 4 He and 3 He, for which the reasonable values 19 MeV 
( 4 He), 3.8 MeV ( 3 He) have been obtained without any adjustable parameters 

III. RESULTS 

A systematic dibaryon search has been done in the u, d and s three flavor world. Both channel coupling and flavor 
symmetry breaking effects have been taken into account. The interesting dibaryon states are listed in Tables 2a, b. 
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Table 2a lists 'Model A' results, where an SU(3) flavor symmetric WF is used and the one body energy difference due 
to m s ^ m have also been neglected. The two body matrix elements are assumed to obey the (empirically consistent) 
relations 

(sl\H l3 \ S l) = -(UlHylll) 

(aa\Hij\aa) = ^(U\H l3 \ll) (17) 

(ll\Hij\ll) = A is calculated using the massless quark WF (Eq.6) operator matrix element multiplied by the overall 
constant such that g 2 s {m — 0)A= —1.350 GcV. The notation sl(ss) means strange-nonstrange (strange-strange) two 
body matrix elements. Table 2b lists the 'Model B' results where the one-body energy for the u,d and s quarks and 
the two-body matrix elements are both calculated with the WF obtained from Eq.(5) with m=0 and m s = 307 MeV, 
respectively. The overall constants, g 2 (0)A= —1.350 GeV and g s (0)g s (m s )A and g 2 (m s )A appear in (ll\Hij\ll), 
(sl\Hij\sl) and (ss\Hij\ss) respectively to obtain the full flavor dependence of the two body matrix elements. The 
value g s (m s ) = <? s (0)/1.26 determined by the ground state baryon masses has been used. 

In principle we should use an SU(3) flavor symmetry breaking WF to study the flavor symmetry breaking effects. 
We use these two model results to show the uncertainty of the flavor symmetry breaking effect in our model. The 
largest difference is for the a — (—600) state where the model A and B results have 140 MeV difference. Except for 
this extreme strangeness case, the largest difference is about 60 MeV corresponding to a = (—3^2). The channel 
coupling effects are similar to the nonrelativistic case. 

The most prominent feature is that the relativistic model yields dibaryon masses in the whole u, d, s three flavor 
world quite similar to those of the nonrelativistic model. This is shown in Table 3, where NA, NB, RA and RB stand 
for nonrelativistic color screening version A and B |6| and the relativistic model A and B (with channel coupling and 
flavor symmetry breaking, ccb) results, respectively. The similarity is true also for states omitted from these tables. 
This might be considered surprising as the relativistic and nonrelativistic model are not only different in kinematics 
but also in the details of the confinement mechanism. We take this similarity of results as a mutual confirmation of 
the stability to the model details of the estimated dibaryon masses in these two models. 

Comparing the relativistic (R) and nonrelativistic (N) model estimates on the dibaryon masses more carefully, one 
finds there are differences fluctuating from case to case. 

In the nonstrange sector, the R model masses are higher than the N ones. In the R model, the deuteron is unbound 
and the model a = (003) state is ~60 MeV higher than the N one. For the strangeness -1 sector, the N and R models 
give almost the same dibaryon masses. 

Beginning at strangeness -1 and generally increasing with increasing strangeness, the R model gives larger binding. 
Except for the a — (-220) state, where the R model mass is 10 MeV higher than the N model mass, in all the other 
cases, R masses are smaller than N masses, ranging from -25 MeV for the a — (-400) state to -220 MeV for the a 
= (-600) and a — (-3^2) states. The a — (-600) state is a special example, where the two versions of the R model 
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themselves have a 140 MeV difference, while all others differ by less than ±60 MeV. However it is interesting to note 
that skyrmion model practitioners have predicted the same strangeness -6 dibaryon candidate jl4j . 

This general trend of differences between the R and N models may be understood in the following way: In the 
N model fit to the octet and decuplet baryon masses, the theoretical octet masses are larger than experiment while 
the decuplet masses are smaller. This indicates that the color magnetic interaction decreases (theoretically) more 
than warranted from the nonstrange to the strange quark case, and is due to the relatively large strange quark mass 
needed for the overall fit to the data. In addition, the R model calculations for dibaryon masses use variational upper 
bounds for the kinetic energies, and so should lead to overestimates in all cases. Finally, the two versions of the R 
model differ from each other more with increasing strangeness since they treat the flavor symmetry breaking in the 
one-body matrix elements differently, as described earlier. 

Excluding a few special cases, the mass differences, among all of these model estimates are in the neighborhood of 
±50 MeV, on average. We take this as the uncertainty of the present dibaryon mass estimate, using the R and N 
models altogether. 

Although our primary interest has been in the a=(003), the R model predicts an H particle state very similar 
to that of Jaffe |ffl], which has already been extensively searched for. Experimental results jl5| on doubly strange 
hypernuclei cast doubt on its existence as a bound state, (although the issue of the relative binding energy of an H in 
a nucleus vs that of two A's has not, to our knowledge, yet been addressed seriously). In the R model, it is ~60 MeV 
lower than the AA threshold, and its wave function is almost a pure flavor singlet as is Jaffe's. However, since the N 
model finds the H to be unbound by 35 MeV, this state is clearly sensitive to dynamical details. Therefore, in our 
models, it cannot definitively theoretically be concluded that the H exists as a bound state. 

A preliminary relativistic calculation of the a=(003) state was reported in Ref. where this state is denoted d*. 
This AA state has the largest binding energy in the it, d and s three flavor world. Because it is a nonstrange state, 
the mass estimated by our model (both relativistic and nonrelativistic) should be more reliable. Even if we take the 
highest estimated mass, it is still 15 MeV lower than NN7T7T threshold. Therefore the d* is quite possibly a narrow 
dibaryon resonance. For comparison, a small-hard-core radius, large NNp-coupling meson exchange model obtained 
a similar state Jl6|] . 

The strange sector also has a few dibaryon candidates. We note especially the spin 3 states, which all have a larger 
predicted binding energy than the d*. They have been discussed in Ref. we refer the reader to that discussion. 

As to the question of prospects for experimental searches: There have, of course, been many efforts to search for 
Jaffe's H. For the d* state, a preliminary estimate of the 7r ± c? — > 7r ± d* production cross section is of order 0.1 fib |fL3| . 
This may be too small to be detected within the large general pion scattering background. A similar reaction |p"7j| 

7r~ ± 3 He -> n + d* -> n + n + n + ir + (18) 

and its charge conjugate reaction 
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n + + 6 H->p + <T->p + p + p + ir . (19) 

may be more favorable for the detection of d* production. Here one can suppress the background by measuring the 
'spectator' n(p) and the emitted tt + (tt~) in coincidence. Unfortunately, pion beams of the required energy range and 
intensity do not seem to be currently available at any accelerator facility. 

On the one hand, coupling to the 1 = 0, 3 D^ NN channel is expected to be small, obviating a search by this 



method. However on the other hand, an analysis by Lomon |18| suggests that the channel coupling may not be small 
enough to avoid conflict with results from phase shift analyses of scattering data already available, even though the 
experimental points are widely spread in energy. 

Finally, we note that a recent calculation by Wong [|l9| of a proton induced reaction similar to the pion induced one 

p+ 2 H^p + d*^p + p + p + n- (20) 

may produce sufficient signal for observation, again using detection of the 7r~ to suppress backgrounds. Such an 

experiment has been proposed at TR1UMF (Canada) by S. Yen (Spokesman for E772). 
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Table 2a. Model A results: 
An SU(3) flavor symmetric WF is used and the one body 
energy difference due to m s / m has been neglected. 
s is the value of x (in fm.) at which the minimum 
energy of the dibaryon state occurs. 



SI J 




M Q 


V a e s 


Threshold 


001 


scs 


1902 


-8 0.2 1.4 


1878(NN) 




CCS 


1898 


-8 0.2 1.5 




010 


scs 




0.0 3.0 


1878(NN) 




CCS 




0.0 3.0 




003 


scs 


2143 


-349 1.0 1.3 


2464(AA) 










2158(NN7T7r) 


1 i o 

-1±3 


scb 


2309 


-334 1.0 1.3 


2617(A£*) 




ccb 


2309 


-334 1.0 1.3 2335(NA7T7r) 


-1|0 


scb 


2145 


-9 0.2 1.6 


2132(NS) 




ccb 


2145 


-9 0.2 1.6 




-200 


scb 


2110 


-230 0.4 0.8 


2231 (AA) 




ccb 


2173 


-242 0.3 0.7 


O AO 

-zOz 


scb 


007A 


-235 0.6 1.1 


2472(NS*) 




CCD 


2371 


-240 0.5 1.0 


2397(NStt) 


-zlo 


scb 




-322 1.0 1.3 


2765(AS*) 




ccb 


Z^± ( u 


-317 1.0 1.3 


2690(AStt) 


-220 


scb 


2396 


-12 0.2 1.5 


2386(SS) 




ccb 


2396 


-12 0.2 1.5 




-3§3 


scb 


2616 


-312 1.0 1.3 


2904(AO) 




ccb 


2634 


-302 1.0 1.3 2788(AS*tt) 


-3|1 


scb 


2494 


-49 0.4 1.2 


2511(ES) 




ccb 


2496 


-50 0.4 1.2 




-3±2 


scb 


2405 


-255 0.6 1.0 


2611(Nfi) 




ccb 


2481 


-252 0.4 0.9 


2574(AStt) 


-3±1 


scb 


2350 


-132 0.6 1.0 


2434(AS) 




ccb 


2456 


-139 0.4 0.9 




-400 


scb 


2622 


-33 0.3 1.5 


2636(SS) 




ccb 


2623 


-33 0.3 1.5 




-600 


scb 


3232 


-131 0.5 1.4 


3345(017) 



scs - single channel only, flavor symmetry 
ccs - with channel coupling, flavor symmetry 
scb - single channel only, broken flavor symmetry 
ccb - with channel coupling, broken flavor symmetry 
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Table 2b. Model B results: 
One-body energy and the two-body matrix elements are both 
calculated with the WF solutions for m=0 and m s — 307 MeV. 



SIJ 




M Q 


V a e s 


Threshold 


001 


scs 


1902 


-8 0.2 1.4 


1878(NN) 




CCS 


1898 


-8 0.2 1.5 




010 


scs 




0.0 3.0 


1878(NN) 




CCS 





0.0 3.0 




003 


scs 


2143 


-349 1.0 1.3 


2464(AA) 










2158(NN7T7r) 


-1±3 

2 


scb 


2293 


-354 0.9 1.2 


2617(A£*) 




ccb 


2293 


-354 0.9 1.2 2335(NA7T7r) 


-l|0 

2 


scb 


2145 


-9 0.1 1.6 


2132(NS) 




ccb 


2145 


-9 0.1 1.6 




-200 


scb 


2071 


-269 0.4 0.8 


2231 (AA) 




ccb 


2171 


-274 0.2 0.6 




-202 


scb 


2246 


-277 0.5 1.0 


2472(NS*) 




ccb 


2332 


-286 0.4 0.9 


2397(NStt) 


-213 


scb 


2432 


-363 0.9 1.2 


2765(AS*) 




ccb 


2439 


-358 1.0 1.2 


2690(AStt) 


-220 


scb 


2393 


-22 0.3 1.3 


2386(SS) 




ccb 


2393 


-22 0.3 1.3 




-3|3 


scb 


2564 


-374 0.8 1.1 


2904(AO) 




ccb 


2581 


-363 0.7 1.1 2788(AS*tt) 


-3|1 


scb 


2459 


-98 0.3 1.0 


2511(SS) 




ccb 


2462 


-99 0.3 1.0 




-3±2 


scb 


2346 


-326 0.6 0.9 


2611(NO) 




ccb 


2420 


-327 0.4 0.8 


2574(AStt) 


-3±1 


scb 


2296 


-197 0.6 0.9 


2434(AS) 




ccb 


2420 


-210 0.3 0.8 




-400 


scb 


2591 


-71 0.3 1.3 


2636(SS) 




ccb 


2591 


-71 0.3 1.3 




-600 


scb 


3093 


-281 0.5 1.1 


3345(017) 



Second column labels are same as in Table 2A. 
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Table 3. The Comparison of the masses of six-quark system between different versions of model. 



SIJ 


\T A 1 n £ — 'Z 

NA: v = 1.6/m 


1\TT~) C\ O £ — 'A 

NB: v = 0.6,/m 


RA 


RB 




M Q e s 


M a e s 


M Q e s 


M Q e s 


001 
010 
003 


1874 0.2 1.3 
1885 0.1 1.6 
2084 1.0 1.2 


1879 0.3 1.3 
1890 0.2 1.4 
2073 1.0 1.2 


1898 0.2 1.5 
0.0 3.0 
2143 1.0 1.3 


1898 0.2 1.5 
0.0 3.0 
2143 1.0 1.3 


-1±3 
-lfO 


2292 1.0 1.1 
2131 0.2 1.4 


2285 1.0 1.1 
2127 0.2 1.4 


2309 1.0 1.3 
2145 0.2 1.6 


2293 0.9 1.2 
2145 0.1 1.6 


-200 
-202 
-216 
-zzU 


2255 1.0 0.6 

2456 1.0 0.8 

2508 1.0 1.0 
oqq^ no 1 q 


2257 1.0 0.6 

2455 1.0 0.8 

z495 l.U 1.1 
oqco no 1 Q 


2173 0.3 0.7 

fill f-1 /\ r~ 1 (\ 

2371 0.5 1.0 

O A '~7f i 1 A 1 O 

2476 1.0 1.3 
zoyo u.z i.o 


2171 0.2 0.6 
2332 0.4 0.9 

O/ion 1 n 1 o 

24o9 l.U 1.2 
zoyo U.o l.o 


-3|3 
-3§ 1 
-3±2 
-3±1 


2711 1.0 0.9 
2509 0.4 1.0 
2669 1.0 0.6 
2560 1.0 0.6 


2697 1.0 1.0 
2506 0.4 1.0 
2670 1.0 0.6 
2562 1.0 0.6 


2634 1.0 1.3 
2496 0.4 1.2 
2481 0.4 0.9 
2456 0.4 0.9 


2581 0.7 1.1 
2462 0.3 1.0 
2420 0.4 0.8 
2420 0.3 0.8 


-400 


2632 0.2 1.1 


2630 0.3 1.1 


2623 0.3 1.5 


2591 0.3 1.3 


-600 


3375 1.0 0.6 


3376 1.0 0.6 


3232 0.5 1.4 


3093 0.5 1.1 
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